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The problem of deriving from microscopic theory a Ginzburg-Landau free energy functional to 
describe the Peierls or charge-density-wave transition in quasi-one-dimensional materials is consid- 
ered. Particular attention is given to how the thermal lattice motion affects the electronic states. 
Near the transition temperature the thermal lattice motion produces a pseudogap in the density of 
states at the Fermi level. Perturbation theory diverges and the traditional quasi-particle or Fermi 
• liquid picture breaks down. The pseudogap causes a significant modification of the coefficients in 

\ the Ginzburg-Landau functional from their values in the rigid lattice approximation, which neglects 

^\ ■ the effect of the thermal lattice motion. 
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I. INTRODUCTION 
A. Motivation 



A wide range of quasi-one-dimensional materials undergo a structural transition, known as the Peierls or charge- 
density-wave (CDW) transition, as the temperature is loweredETQ. A periodic lattice distortion, with wave vector, 
2kp, twice that of the Fermi wavevector, develops along the chains. Anomalies are seen in the electronic properties 
due to the opening of an energy gap over the Fermi surface. 

Over the past decade, due to the development of high-quality samples and higher resolution experimental techniques, 
new data has become available which allows a quantitative comparison of experiment with theory. The most widely 
studied material is the blue bronze, K0.3M0O3. There is a well-defined three-dimensional transition-at Tp = 183 K and 

■ careful measurements have been made of thermodynamic anomaliesB and CDW coherence lengthsa at the transition. 
The critical region, estimated from the Ginzburg criterion!!] is only a few percent of the transition temperature and 

O ■ so the transition should be described by an anisotropic three-dimensional Ginzburg-Landau free energy functional, 
' except close to the transition temperature. The challenge is to derive from a microscopic theory the coefficients in the 

^ , Ginzburg-Landau free energy so a quantitative comparison can be made between theory and experiment. Inspiration 
is provided by the case of superconductivity. The superconducting transition is well described by Ginzburg-Landau 
theory and the coefficients can be calculated from BCS theory! and depend on microscopic parameters such as the 

■ normal state density of states, Debye frequency, and the electron-phonon coupling. This program is so successful that 
one can even consider refinements to BCS theory, such as strong coupling effects, in order to get better agreement 
between experiment and theoryEI. However, the problem of the CDW transition is more difficult because of the large 
fluctuations due to the quasi-one dimensonality. 



B. Ginzburg-Landau theory 

The Peierls transition is described by an order parameter which is proportional to the 2kp lattice distortion along 
the chains. The order parameter is complex if the lattice distortion is incommensurate with the lattice. For a 
commensurate lattice distortion (e.g., a half- filled band) the order parameter is real. I recently considered the general 
problem of Ginzburg-Landau theory fpc a three-dimensional phase transition, described by a complex order parameter, 
in a system of weakly coupled chainai3. The key results of that study are now summarized, partly to put this paper 
in a broader context. 

The Ginzburg-Landau free energy functional F\[<j>\ for a single chain with a complex order parameter <f>(z), where 
z is the co-ordinate along the chain, is 
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Near the single chain mean-field transition temperature To the second-order coefficient a(T) can be written 

T 



a(T) = a! \ I 
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Due to fluctuations in the order parameter this one-dimensional system cannot develop long-range order at finite 
temperatureO'E3. To describe a finite-temperature phase transition, consider a set of weakly interacting chains. If 
4>i(z) is the order parameter on the i-th chain the free energy functional for the system is 



(3) 



<i,j> ' 



where J describes the interchain interactions between nearest neighbours. A mean-field treatment of this functional 
will only give accurate results if the width of the three-dimensional critical region is much smaller than T - This 
requires that the width of the one-dimensional critical region Atm = (Wo) 2 / 3 /a'c 1 / 3 be sufficiently small that 



At 1£ > < 



J 



2/3 



(4) 



If this is not the case one can integrate out the one-dimensional fluctuations to derive a new Ginzburg-Landau 
functional with rcnormalized coefficients, 
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where a x and a y are the lattice constants perpendicular to the chains. The new order parameter y, z), is 
proportional to the average of 4>i(z) over neighbouring chains. The three-dimensional mean- field temperature T^£> is 
defined as the temperature at which the the coefficient A(T) changes sign. Close to T^u 



A = A 
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The transition temperature T30 and the coefficients A' , B, C x , C y , and C z can be written in terms of the interchain 
interaction J and the coefficients a, b, and c of a single chain. The coefficients in (|^) determine measurable quantities 
associated with the transition such as the specific heat jump, coherence lengths and width of the critical region. 

Most of the physics is determined by a single dimensionless parameter 
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which is a measure of the fluctuations along a single chain. It was assumed that the coefficients a, b, and c were 
independent of temperature and the measurable quantities at the transition were determined as a function of the 
interchain coupling. The transition temperature increases as the interchain coupling increases. The coherence length 
and specific heat jump depends only on the single chain coherence length, £ = (c/laQ 1 / 2 , and the interchain coupling. 
The width of the critical region, estimated from the Ginzburg criterion, was virtually parameter independent, being 
about 5-8 per cent of the transition temperature for a tetragonal crystal. Such a narrow critical region is consistent 
with experiment, and shows that Ginzburg-Landau theory should be valid over a broad temperature range. 

This paper uses a simple model to demonstrate some of the difficulties involved in deriving the coefficients a, b, c, 
and J from a realistic microscopic theory. 



C. Microscopic theory 

The basic physics of quasi-one-dimensional CDW materials is believed to be described by a Hamiltonian due to 
FrohlichuJ which describes electrons with a linear coupling to phonons. Even in one dimension this isr-a [highly 
non-trivial many-body system and must treated by some approximation scheme. The simplest treatment£3ii3 is a 
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rigid-lattice one in which the phonons associated with the lattice distortion are treated in the mean-field approximation 
and the zerp-point and thermal lattice motions are neglected. The resulting theory is mathematically identical to 
BCS theoryta. An energy gap opens at the Fermi surface at a temperature Trl ~ l.lAEpe where Ep is the 
Fermi energy and A is the dimensionless electron-phonon coupling. Trl is related to the zero-temperature energy gap 
A„,.(0) by 



A RL (0) = 1.76k B T RL . 

In this approximation the coefficients in the single-chain Ginzburg - Landau free energy functional (IT]) areEi 



(irl(T) 



1 m(J-) 

Trl J 



7TV p 



(8) 



(9) 



b RL (T) 



1 7C(3) 

TTVp (AnT) 2 



(10) 



c RL {T) 



1 7C(3) V 



TTVp (AnT) 2 



where vp is the Fermi velocity and £(3) is the Riemann zeta fun 
to the chains (see (pTl)) then the interchain coupling is given byH 



(11) 



If 4t± is the electronic bandwidth perpendicular 



Jrl = 
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It might be hoped that the transition in real materials can be described by the mean-field theory of the functional 
(|lj) with the coefficients (^-^JJ). However, this is not the case for several reasons_£i) The width of the critical regime 
given by the one-dimensional Ginzburg criterionE£l is very large: AtiD = 0-8l3, suggesting that fluctuations are 
important because condition (^) is not satisfied, (ii) A rigid-lattice, treatment predicts a HietsJlic density of states fU, all 
temperatures above Trl- In contrast, magnetic susceptibilitylijiHJ, optical conductivityEftll, and photoemissionE§~E2l 
measurements suggest that there is a gap or pseudogap in the density of states for a broad temperature range above 
Tp. (iii) The transition temperature, specific heat jump, and coherence lengths are inconsistent with rigid lattice 
predictions (Table |). This failure should not be surprising given that recent work has shown that in the three- 
dimensionally ordered Peierls state the zero-point .and .lijjexmal lattice motions must be taken into account to obtain 
a quantitative description of the optical propertiesScj'EJTJ. 

The next level of approximation is to use the coefficients (|^|l0|) and take into account the intrachain order parameter 
fluctuations and the interphain-coupling and use results similar to those in ReferencesEl This is the approach that 
has been taken previouslyc3'OE3. There are two problems with this approach. First, if the dimensionless parameter 
k, given by (0), is evaluated using the expressions (§|ll]) the result is 
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Hence, tri p .. temp erature dependence is quite different from the dependence k ~ T 2 that was assumed in 
Referencesli^O^j and the analysis there needs to be modified. The second and more serious problem is one of 
self-consistency. The coefficients a, 6, and c are calculated neglecting fluctuations in the order parameter which will 
modify the electronic properties which in turn will modify the coefficients. In this paper a simple model is used to 
demonstrate that the fluctuations have a significant-effect on the single-chain coefficients. 

An alternative microscopic theory, due to Schulzt3, and which takes into account fluctuations in only the phase of 
the order parameter is briefly reviewed in Appendix |A|. 



D. Overview 



Discrepancies between phonon rigid-lattice theoriL.and the observed properties of the Peierls state well below the 
transition temperature Tp were recently resolvecBES by taking into account the effect of the zero-point and thermal 
lattice motion on the electronic properties. It was shown that the lattice fluctuations have an effect similar to a 
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Gaussian random potential. This mapping breaks down near the transition temperature because of the phonon 
dispersion due to the softening of the phonons near 2kp . In this paper this dispersion is taken into account and the 
effect of the large thermal lattice motion near the transition temperature is studied. 

The thermal lattice motion has the same effect on the electronic properties as a static random potential with finite 
correlation length. Close to the transition temperature, the problem reduces to a simple model, corresponding_to a 
single classical phonon, which can be treated exactly (Section This model was first studied by SadovskiO. It 
was recently used in the description of the destruction of spin-density-wave states by high magnetic fieldsEa. The 
one-electron Green's function is calculated in Section II] . There is a pseudogap in the density of states (Fig [j]) . The 
complexity of this simple model is indicated by two non-trivial many-body effects: (i) Perturbation theory diverges but 
is Borel summable.|-(jii) The traditional quasi-particle picture breaks down (Figure ||), reminiscent of behaviour seen 
in Luttinger liquidsa. To illustrate that calculations based on perturbation theory can be unreliable it is shown that, 
a predicted scaling relation between the specific heat and the temperature derivative of the magnetic susceptibility^ 
does not hold if the exact, rather than approximate, density of states is used in the calculation. Using this model 
the coefficients a, b, and c are calculated in Appendix [§[ The coefficients deviate significantly from the rigid-lattice 
values if the pseudogap is comparable to or larger than the transition temperature. In Section VI experimental data 
is used to estimate the pseudogap in K0.3M0O3. 



E. Previous work on fluctuations and the pseudogap 



To put this paper in context some important earlier work is briefly reviewed. 

Lee, Rice and AndersonEfl considered how fluctuations in the order parameter produce a pseudogap in the density 
of states. It is important to be aware of the assumptions made in their calculation. Although their results describe 
much of the physics on a qualitative level, for the reasons described below, their results cannot be expected to 
give a quantitative description of the density of states near the CDW transition. The starting point of Lee, Rice, 
and Anderson was the one-dimensional Ginzburg-Landau functional (|l|) with a real order parameter a»d with the 
coefficients derived from rigid-lattice theory (see equations (|s|-|lT|) ) . Earlier, Scalapino, Sears, and Ferrel£3 evaluated 
the correlation length £y (T) for one-dimensional Ginzburg-Landau theory with an exact treatment of the fluctuations 
in the order parameter; £11 (T) only diverges as T — > 0. The results of this calculation were used by Lee, Rice, and 
Anderson as input in a random potential with correlations given by 

< A(z)A(z') >= A RL (Tf exp(- \z-z\ /£„ (T)) (14) 

where A b.l (T) is the rigid-lattice (BCS) order parameter and the average is over the thermal fluctuations of the order 
parameter. The electronic Green's function was calculated using equation ( |l4| ) and a formula originally used for liquid 
metals (essentially, second-order perturbation theory for the random potential). They found a gradual appearance 
of a gap as the temperature decreased. For Tp < Q.25Trl, an absolute gap of magnitude A#l(0) appears. Lee, 
Rice and Anderson suggested that a three-dimensonal transition occurs for Tp ~ 0.25Tpl based on the temperature 
at which £11 (T) becomes extremely large. There are several problems with trying to use these results to give a 
quantitative description of the CDW transtion, because of the following assumptions, (i) A real order parameter. 
Most CDW transitions are described by a complex order parameter, for which quantitatively distinct behaviour 
occurs. For example, the transition to very lacge correlation lengths for Tp ~ 0.25Trl does not occur for a complex 
order parameter. (See Figure 6 in ReferencetS) . (ii) Rigid- lattice coefficients. It is shown in this paper that the 
pseudogap due to the thermal lattice motion causes the Ginzburg-Landau coefficients to deviate significantly from 
their rigid-lattice values (Figure 3). (iii) Perturbation theory. It is demonstrated in this paper that this is unreliable. 
In particular as £\\{T) — > 00 in (14) only a pseudogap rather than an absolute gap develops in the density of states 
(Figure|). „ 

Rice and StrassleiH2l calculated the contribution of the phonon fluctuations to the electronic self energy in the Migdal 
approximation, i.e., second-order perturbation theory. Interchain interactions were included through an anisotropic 
phonon dispersion. They found a pseudogap in the density of states above the transition temperature. At Tp there 
is an absolute gap whose magnitude is determined by the electron-phonon coupling and the interchain interactions. 
They equated the observed transition temperature with the single-chain mean-field transition temperature To which 
they found to be significantly reduced below the rigid-lattice value Trl and to vanish as the interchain coupling 
vanishes. 

In the limit .of weak interchain interactions the analytic form of the density of states is identical to that of Lee, Rice, 
and AndersonEj. However, it is not commonly appreciated that the origin of the pseudogap in the two calculations 
is quite different. The magnitude of the Rice and Strassler pseudogap is proportional to the thermal lattice motion 



(compare Section II A), while the pseudogap studied by Lee, Rice, and Anderson pseudogap is by assumption equal 
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to the rigid-lattice gap Arl(T). Calculations similar tq-that of Rice and-Strassler have been performed by Bjelis and 
BarisicCJ, Suzumura and KuriharaCJ, Patton and Shamcfl, and ChandraEI The main problem with these calculations 
are that they are based on perturbation theory. 

II. MODEL HAMILTONIAN 

The starting point for this paper is the following one-dimensional model. The states in an electron gas with Fermi 
velocity vp are described by spinors ^(z). The upper and lower components describe left and right moving electrons, 
respectively. The phonons are described by the field 



AW = *E ^ + (15) 

where b s destroys a phonon of momentum s and frequency lu s and g is the linear electron-phonon coupling. The 
dimensionless electron-phonon coupling A is defined by 

A = 2g 2 a z /nvpu}Q (16) 

where a z is the lattice constant along the chains. The electronic part of the Hamiltonian isEl 



H el = J dz¥(z) 



Q 1 

iv F a 3 — + -(A(z)<t+ + A(z)V_) 
az 2 



*(*) (17) 



where 03 and o± = a\ ± ioi are Pauli matrices. 

This paper focuses on the following model where A(z) is replaced with a random potential with zero mean and 
finite length correlations 

<A(z))=0 <A(z)A(z7)^ 2 exp(-|z-z'|/e||)- (18) 

£|| is the CDW correlation length along the chains. In most of this paper ip will be treated as a parameter. It is 
central to this paper, being a measure of the thermal lattice motion and a measure of the pseudogap in the density 
of states. This paper focuses on behaviour near Tp and so the limit ^hi/j/vf — » oo is taken. A rough argument is now 
given to justify using this model to describe thermal lattice motion near the phase transition. 

A. Thermal lattice motion 

In rigid-lattice theory A(z) is replaced by its expectation value (A(z)) = A D . Tq rgo-beyond this the effect of the 
quantum and thermal lattice fluctuations in the Peierls state was recently modelledLJOcj by treating A(z) as a static 
random potential with mean A and correlations 

(A(z)A(z')*) = A3 + 7<5(z-z') (19) 

where 

7 = -TrXvpUikp coth (^f ) • (20) 

This model is expected to be reliable except near the transition temperature where there is significant dispersion in 
the phonons. This dispersion is now taken into account. 

Near the transition temperature the phonons can be treated classically since in most materials the frequencies of 
the phonons-adth wavevector near 2kp are much smaller than the transition temperature (Table |TJ). Following Rice 
and StasslerEll renormalized phonon frequencies fi(q, T) are used in the expression for the correlations of the random 
potential 

(AWAMVA^Ej^^. (21) 
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At the level of the Gaussian approximation the phonon dispersion relation can be written in the form 

Q(q, Tf = n(T) 2 (1 + (q - 2k F ) 2 ^ (T) 2 ) . (22) 
Evaluating ( |2l| ) then gives ([l8| ) where 

Note that this expression togehter with is then quite different from (O) used by Lee, Rice, and AndersoriEI. In 
the limit ^ — ► 0, i.e., the phonons become dispersionless and the sum in (|2l|) becomes a delta function and giving 
(H) (with A = 0) and ©. 

The rms fluctuations Su in the positions of the atoms due to thermal lattice motion is related to the Debye- Waller 
factor and given by 

(Su) 2 = fcrV i— - (24) 

This is related to i/i — {2MuiQ) 1 / 2 g5u. Hence ip is proportional to the thermal lattice motion. 
If ip is defined by (|23|) it diverges as T — > Tp because 



as T-> T P , fi(T) -> 0, ^||(T) -> oo with n(T)^(T) finite. (25) 
However, in a real crystal the phonons are three-dimensional and the thermal lattice motion is finite. Define 

^ = A.T^£-^ (26) 

and write the three-dimensional dispersion relation (for a tetragonal crystal) in the form 

n(q,T) 2 = n(T) 2 (l + («„ - 2fc F ) 2 e,|(T) 2 + (q ± - Q ± ) 2 ^(T) 2 ) (27) 

where Q = (Q±, 2fcp) is the nesting vector associated with the three-dimensional CDW transition (see equation (j56|)). 
Due to the quasi-one-dimensionality of the crystal the dispersion perpendicular to the chains is small and £j_ <C 
Let (hn denote the lattice constant perpendicular to the chains. Performing the integral over the wave vector in ( |2q ) 



givescJ 

^ = A,T^°V alvF 



) it mT m2 \Vl + (Pc^(T)) 2 -l] (28) 



\Q(T)J n%(T)£ ± (Ty 

where p c is a wavevector cutoff perpendicular to the chains. If p c = Tr/a x this expression reduces to ( |23| ) in the 
one-dimensional limit £± <C a x . Near the transition, £_i_(T) — > oo, giving 

From ( |25| ) and the fact that £|| (T)/£j_(T) is finite it follows that ip is finite as T — > Tp. Note that the magnitude of 
this quantity is dependent on the choice of the momentum cutoff p c . The above treatment is quite similar to Schulz's 
discussion of fluctuations in the order parameter in the Gaussian approximationlia 

Although the expressions ((2^) , ( |28| ) , and (^9| ) for if) in the different regimes look very different ip is actually weakly 
temperature dependent and does not vary much in magnitude. To see this ( p8| ) can be written as 

1 



^ = (<A(Tp)) 2 — — - ^l + ( Pc U(T)) 2 -l (30) 
PcLl(T) L J 

The postfactor is a slowly varying function of p c £±(T). Since well above Tp, p c £±(T) ~ 1 (e.g., for K0.3M0O3 
£_i_(300K) ~ the postfactor does not vary by more than a factor of two although /9 c £j_(T) varies by several 

orders of magnitude. Johnston et al.E2l used a crude method of estimating the pseudogap and found it to be weakly 
temperature dependent above Tp for K0.3M0O3. 



6 



B. Solution of the model 



Sadovskiio solved the one-dimensional model (|T^) and ( |l8| ) exactly. He calculated the one-electron. Green's function 
in terms of a continued fraction by finding a recursion relation satisfied by the self energy. He foundEZI that the Green's 
function reduced to a simple analytic form in the limit of large correlation lengths (£|| 3> vf /ip)- This can be seen by 
the following rough argument. In the limit £|| — > oo the moments of the random potential A(z) are independent of 
position: 

<A(*)> = <A(z)A(z')*) =^ 2 . (31) 

This means that the random potential has only one non-zero Fourier component, i.e., the one with zero wavevector. 

The potential can be written A(z) = vip where v is a complex random variable with a Gaussian distribution. 
Averages over the random potential can then be written 

(A[A(z)}) = J ^l e - m *A[vi>). (32) 

It is then a straight-forward exercise to evaluate the averages of different electronic Green's functions. 



III. ONE-ELECTRON GREEN'S FUNCTION NEAR T P 



The matrix Matsubara Green's function, defined at the Matsubara energies e n = (2n + l)irT, for the Hamiltonian 
(0) with (|l]) is 



G(ie„,k) 



dvdv* 



*G(ie n ,k,v) 



where 



G(ie n ,k,v) 



— (ie n — kvF&3 — ip(va+ + v*a-)) 
+ {kvp) 2 + vv*ip 2 



(33) 



(34) 



is the matrix Green's function for the Hamiltonian ( |17| ) with A(z) = vtp. The off-diagonal (anomalous) terms vanish 
when the integral over v is performed indicating there is no long range order. The integral over the phase of v can be 
performed and variables changed to <p = vv* and obtain 



G(ie n , k) = -{ien - kv F cr 3 ) / dip 



<?n + (kv F ) 2 + ip^ 



(35) 



SadovskiS obtained the same expression by diagrammatic summation. For the case of a half-filled band v is strictly 
real and the resulting expressions are the same as those obtained by Wonneberger and Lautenschlageim Expanding 
( |35| ) in powers of ip gives 



G(ie n ,k) = G (ie n ,k) d(pe~ v ^ 
J o n=0 



{kv F y 



(36) 



where Go = (ie„ — kvpa^) 1 is the free-electron Green's function. Performing the integral over ip gives 



G(ie n , k) = G {ie n , k) : 



(kv f y 



(37) 



This is a divergent series and asymptotic, 
perturbation theory as used in References! 
iHE. 



aion. However, it is Borel summable0. This divergence suggests that 
3 may give unreliable results. This can be seen in Figure |l| and Section 
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A. Density of States 



The electronic density of states is calculated directly from the imaginary part of the one-electron Green's function 
(|35|). The result is 

f°° E 

p{E)= Po dtpe-v— ——9 (| E | 2 -^V 2 ) (38) 

Jo [hi' 1 — ipip z \ ' 

, fE\ 2 ^ „{E 



= ^M-[^) Jerfi^j (39) 

where p D = 1 /nvp is the free-electron density of states and erfi is the error function of imaginary argument^. Figure |l| 
shows the energy dependence of the density of states. It vanishes at zero energy (the Fermi energy) and is suppressed 
over an energy range of order ip, i.e., there is a pseudogap. It has the asymptotic behavior: 

p(E) ~ 2 Po { E -f for£<V (40) 

V 

p(E) ~ p for E > ip 

Figure |l| shows that th|6-exact result ( |39| ) (solid line) deviates significantly from the result of second-order perturbation 
theory in ReferencesE-scil (dashed linejf 

p{E)= PoW ^0{E*-^) (41) 

This latter form has been assumed in much earlier workEZI 

The above expressions for the density of states are all for infinite correlation length (^nip/vp — * oo), i.e., very close 
to the three-dimensional transition temperature Tp. What happens above Tp as the intrachain correlation length 
decreases? This problem was considered in detail by Sadovskfica. (He calculated the density of states for the random 
potential (18) with finite £m exactly). As the correlation length decreases the density of states at the Fermi energy 



increases, i.e., the pseudogap fills in. How quickly this happens depends on the dimensionless ratio vr r/ (i££ i lL (See 
equation below and Figures 5 and 6 in Reference!!) . Sadovskii showed that perturbation theoryOLJ @ only 
gives reliable results for \E\ < -0 when £|| < vf/"4>, i.e., well above Tp. 

What happens below Tp as the intrachain correlation length decreases? In Referencetil it was shown that in the 
three-dimensionally ordered Peierls state, well below Tp, there is an absolute gap with a subgap tail that increases 
substantially as the temperature becomes larger than the phonon frequency. A smooth crossover to the pseudogap 
discussed here is expected. It is an open problem to construct a single theory that can describe the density of states 
over the complete temperature range. 

B. Spectral Function 

The spectral function for right moving electrons of momentum k is given by 
A(k, E) = — Im G n (k, E + irj) 

TT 



dtpe 



-tp 



o 



6 (e — ^{kv F ) 2 +^ 2 ) + 6 (E + ^(kv F ) 2 +<^ 2 ) 



El e.J (kV ^- E2 )e(E>-(k VF r) (42) 



ip \ ip 

where the momentum k is relative to the Fermi momentum kp . Note that this form is very different from the Lorentzian 
form associated with the quasi-particle picture and perturbation theoryEj. The spectral function is asymmetrical, 
very broad, and has a significant high energy tail. Figure ^| shows how the quasi-particle weight is reduced near 
the Fermi momenta, i.e., the quasi-particles are not well defined. This was first pointed out by Wonneberger and 
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LautenschlagerEJ for the corresponding model for a half-filled band. This is strictly a non-perturbative effect. In 
perturbation theory the quasi-particles are well defined. This breakdown of the quasi-particle picture is similar to the 
properties of a Luttinger liquidtS. 

The momentum distribution function n(k) at T = for right moving electrons is given by 

n(k) = / dEA(k, E) = - 

J-oo 2 

where erf is the error function. The inset to Figure^ shows how the momentum distribution n(k) at T — is smeared 
over a momentum range 5k ~ i/i/vp. The absence of a step at k — kp indicates that there is no clearly defined Fermi 
surface. However, this is not like in a Luttinger liquid, but solely due to disorder. In fact, in an ordinary metal with 
mean free path I similar behaviour is seen; disorder smears out n(k) over a momentum range 6k ~ \Ji, 



1 J* (tO exp 




(1 - erf 



(43) 



C. Electronic specific heat 

The electronic specific heat C e (T) is related to the density of states p{E) by 

AC 00 Bf 
Ce(T) = -- J dEE 2 P (E)^L (44 ) 

where f(E) is the Fermi-Dirac distribution function. In the absence of a pseudogap C e {T) = =^-pgT = Cq(T). If the 
expression (|3|) is used for the density of states in the presence of .a, pseudogap then C e (T)/Co(T) only depends on 
ijj/T and is shown in Figure |[ A similar result was recently usedEi to explain the temperature dependence of the 
electronic specific heat near a spin-density-wave phase boundary of the organic conductor (TMTSF^CICU. Note that 
when tp ~ T, C e (T) can be slightly larger than Cq(T) because E 2 ^ has a maximum near E ~ T and for E ~ t/i, 
p(E) is larger than p (Figure |l|) . 



D. Pauli Spin Susceptibility 



The Pauli spin susceptibility x(T) is related to the density of states p(E) by 

X(T) = -V 2 B J™ dEp(E)^L (45) 

where f(E) is the Fermi-Dirac distribution function and p,s is a Bohr magnetonB In the absence of a pseudogap 
x(T) = p-gPo = Xo which is independent of temperature. If the expression (^) is used for the density of states in 
the presence of a pseudogap then xC^O/Xo on ly depends on ip/T and is shown in Figure [| This result will be used 
in Section [yj to provide an estimate of the pseudogap in K0.3M0O3. 



E. Chandra's scaling relation 

The effect o£lhermal lattice fluctuations on the temperature dependence of x(^) was ^ IS ^ considered by Lee, Rice, 
and Andersora. They argued that as the temperature is lowered towards Tp the intrachain correlation length in- 
creases, moEe-pf a pseudogap opens in the density of states and x(T) decreases. This problem was recently reconsidered 
by ChandraEl who derived a scaling relation between the derivative dx/dT and the specific heat Cp in the critical 
region. I now repeat the essential features of her argument. She calculated the electronic self energy in the Born 
approximation, taking into account the interchain interactions and the finite mean free path of the electrons. She 



assumed that the pseudogap is much larger than the transition temperature (ip 3> Tp; it will be shown in Section VI 
that this is poor approximation for K0.3M0O3) so that x(T) ~ p%p(0). Chandra also assumed that the temperature 
dependence of the density of states at the Fermi energy is determined solely by the temperature dependence of £|| (T). 
Moreover, based on the Born approximation, she found 

«»~m (46> 







Defining t = \T — Tp\/Tp, then gives the scaling relation 



dT 



,1/2 



C P 



(47) 



dT C|| (T) dT 

where use has been made of the temperature dependence of £|| (T) and Cp in the Gaussian approximationlll. 

This same scaling relation was suggested earlier by Horn, Herman and SalamonFjj. They claimed to have found 
the critical exponent for d\/dT to be -0.5 for TTF-TCNQ. Kwok, Griiner, and BrownEJ claim to have observed a 
scaling between d(Tx)/dT and Cp within a 30 K region about Tp = 183 K for K^MoOs. Howejaer, Mozurkevich has 
argued that the Gaussian approximation is not valid in this temperature rangeEZI. Chung et a^Hrfound that d\/dT 
was comparable to Cp when a background contribution was subtracted from the latter. Brill et a/.El found that \ was 
proportional to the entropy (evaluated from integrating the specific heat) between 140 and 220 K. (This is equivalent 
to a scaling between d\/dT and Cp). They show that this is what is expected if \ an d Cp are derived from a free 
energy functional in which the complete magnetic field dependence is contained in the field dependence of Tp. 

Chandra's derivation of the scaling relationship ( |47| ) is not valid. It depends on ( flrf ) which is a direct result of the 
perturbative treatment of the lattice fluctuations. The exact Green's function calculated by Sadovskiica gives different 
results. He found that for £||(T) ^> vp/ip 



(0.54±0.01)( 



v F 



1^1 CO 



1/2 



(see Figure 6 in Referenceca) rather than (fi"6|). This will give 

d X(T) 3/4 

dT ~ l 



(48) 



(49) 



and so the scaling relation (47) does not hold. It should be stressed that this result assumes ip 3> Tp, a condition 
that is poorly satisfied in most materials (Section [y|) . 



IV. PROPERTIES OF THE GINZBURG-LANDAU COEFFICIENTS 



In Appendix [b| the coefficients a, b, and c in the Ginzburg-Landau free energy (Q) describing the Peierls transition 
are evaluated in the presence of the random potential (31) which is used here to model the thermal lattice motion. 
The calculation is hased on a linked cluster expansion similar to that used to derive the Ginzburg-Landau functional 
for superconductorsE^I. The results are: 



a(T) = — 

TTVf 



In 



Trl) 



d(pe 



(el + (pip* 



^3/2 



(50) 



dipe~ 



\5/2 



7/ 2 



(51) 



(rp\ l 'FT X - 2 



dipe 



(52) 



The integrals over <p in the above expressions can be written in terms of error functions and incomplete gamma 
functionsLil. However, for both numerical and analytical calculations it is actually more convenient to use the expres- 
sions above. As ip —* the above expressions reduce to the rigid-lattice values (§-11). 

Single-chain mean-field transition temperature. To is determined by the temperature at which the second-order 
Ginzburg-Landau coefficient (|5(]) vanishes: 



a(T ) = 0. 



(53) 



This defines relations between Tq/Trl and ip/To, shown in Figure |j (The inset shows To /Trl versus iP/Trl)- The 
pseudogap suppresses the transition temperature. At a crude level, this is because in the presence of a pseudogap 
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opening a gap due to a Peierls distortion causes a smaller decrease in the electronic energy than in the absence of 
a pseudogap. In most materials Tp < OATpx (Table ||) and so the inset of Figure || implies tp ~ Trl which is 
comparable to the zero-temperature gap. Rice and Strasslercll found from second-order perturbation theory that for 
To -C Tm,, tp ~ 1.05Trl. Thus, the single-chain mean-field transition temperature can be quite different from Trl, 
defined by (Q), and often referred to as the mean- field transition temperature, and so no experimental signatures are 
expected at T = Trl. 

Fourth-order coefficient. The ratio of the fourth-order coefficient b to its rigid-lattice value as a function of the ratio 
of the pseudogap tp to the temperature is shown in Figure || Note that b is negative for tp/T > 2.7. This will change 
the nature of the phase transition. One must then include the sixth-order term in the free energy. If it is positive 
(I have calculated it and found it to be positive for this parameter range V4hen the transition will be first order. A 
complete discussion of such a situation is given by Toledano and ToledanoEJ. Imry and Scalapino have discussed the 
effect of one-dimensional fluctuations for this situation!^. At the mean-field level there is a co-existence of phases for 
the temperature range defined by 

0<«Cr)<|g (54) 

where d(T) is the sixth-order coefficient. Hysteresis will be observed in this temperature range. I recently suggested 
that the first-order natupe of the destruction by high magnetic fields of spin-density-wave states in organic conductors 
is due to similar effectsH. If at low temperatures the electron phonon coupling A is varied then tp/Tp^ ~ Ae 1 '*. 
According to the inset of Figure [| there will be a critical coupling below which -the CDW phase will be destroyed. 
This transition will be first order. It is interesting that Altshuler, Ioffe, and MilliscB recently obtained a similar result 
for a two-dimensional Fermi liquid (with a quasi-one-dimensional Fermi surface) using a very different approach. 

However, it should be pointed out that when b is small corrections due to other effects such as a finite correlation 
length and interchain coupling will be important and could make b positive. It is unclear whether this unexpected 
behaviour is only a result of the simplicity of the model or actually is relevant to real materials. The three-dimensional 
transition occurs when the parameter k, defined by (Q), becomes sufficiently smalfcj. Generally this is assumed to 
be due to the temperature becoming sufficiently low. However, I speculate that the transition could alternatively be 
driven by b becoming sufficiently small. The fact that tp ~ (2 — 3)Tp in K0.3M0O3 (Section |V^ ) is consistent with b 
being small. 

The coefficient of the longitudinal gradient term is given by (|52|). It can be shown that c{T)/c RL (T) is a universal 
function of tp/T (see Figure ||) and that the pseudogap reduces the value of c. 

Interchain coupling. Consider a crystal with tetragonal unit cell of dimensions a x x a x x a Zl where the z-axis is 
parallel to the chains. For a tight-binding model the electronic band structure is given by the dispersion relation 

E(k) = —2t±(cos(k x a x ) + cos(k y a x )) — 2t\\ cos(k z a z ). (55) 

Assume the band-structure is highly anisotropic—i.e., t\\ t±. The Fermi velocity vp is defined by vp — 
2t\\a z sin(kpa z ). Horovitz, Gutfreund, and WegeiO have shown that imperfect nesting of the Fermi surface (i.e., 
E(k) ~ —E(k + Q)) occurs for the nesting vector 

Q = (Tr/a x , Tr/a x ,2k F ). (56) 

To calculate the interchain coupling J in the Ginzburg-Landau functional (||) it is assumed that the one-dimensional 
Green's function (|34|) can simply be replaced with the corresponding one with the anisotropic band structure, given 



it ^ 2 



by equation-.(55). The calculation is then essentially identical to the rigid-lattice calculation of Horovitz, Gutfreund, 
and Wegertj and so only the result is given (compare (|l2])): 

■) c(T). (57) 

V F J 

Since the pseudogap reduces the value of the longitudinal coefficient c it will also reduce the interchain coupling. 

V. MEAN-FIELD THEORY OF A SINGLE CHAIN 

The single chain Ginzburg-Landau functional with the coefficients discussed in the previous section is now con- 
sidered. In particular it is shown that the one-dimensional fluctuations can be much smaller than for the functional 
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with the rigid-lattice coefficients. The first step is to consider the temperature dependence of the second-order coef- 
ficient a(T) near To, the mean-field transition temperature. This is difficult because to be realistic the temperature 
dependence of the parameter ip must be includ ed. T his is done at a crude level using the simple model based on the 



discussion of thermal lattice motion in Section If A . This is then used to evaluate a', defined by (|2|), and needed to 
evaluate physical quantities associated with the transition: the specific heat jump, the coherence length, and width 
of the critical region. 

The jump in the specific heat at To is 

(a') 2 

AC 1C . ±J_. (58) 
An important length scale is the coherence length £o, defined by 

1/2 



/ C \ 1 z 

Co = y (59) 



The one-dimensional Ginzburg criterionfl provides an estimate of the temperature range, ATid, over which critical 
fluctuations are important. 

AT 1D ( 6 T ^ 2/3 



A. Self-consistent determination of the pseudogap 

At the level of the Gaussian approximation the phonon dispersion is related to the Ginzburg-Landau coefficients 

by 

n(q,T) 2 = \Lu 2 Q (a(T) + c(T){q-2k F ) 2 + Ja 2 x ( qi _ -Qa_) 2 ). (61) 



Hence the phonon dispersion depends on the pseudogap ip. However, it was shown in Section II A that ip depends on 
the dispersion. Hence, ip must be determined self-consistently. Equation (|2^) gives the dependence of the pseudogap 
at To on the phonon dispersion. Equation ([52]) gives the dependence of the coefficient c(T) on the pseudogap. These 
can be combined with ( |57|) to give 



<kpe 



l=^ 2 E e »/ 73 Tl^,, - (62) 



It follows that ip/T is a universal function of t±/T. 

Dependence of Tq on the interchain interactions. The self-consistent equation for the pseudogap (|6^) can be solved 
simultaneously with the equations for To, and (p7| ) to give the single-chain mean-field transition temperature as 
function of the interchain interactions. The transition temperature is then a-^ponotonic increasing function of the 
interchain hopping. A similar procedure was followed by Rice and Strasslercil. The transition temperature tends 
to zero as the interchain coupling tends to zero—consistent with the fact that there are no finite temperature phase 
transitions in a strictly one-dimensional systemO. 



B. Evaluation of a' 

It is now assumed that the temperature dependence of the pseudogap ip is given implicitly by equation Jo^ ) . Implicit 
differentiation then gives 

d ftp\ ib X(T) , x 

' v ' - V V ' (63) 



dT \T J 2T 2 Y{T) 
where 

dipe~ v 
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Note that since the right-hand side of (|63| ) is positive ip/T is always an increasing function of temperature. A lengthy 
calculation gives 



„' = J_ ( l + U\ T V el r , 2 dipC Z 5/2 ) (66) 



This is large than the rigid-lattice value a' RL = 1/wvf- This enhancement will enhance the specific heat jump (|5 
and reduce the coherence length 



C. Specific heat jump 

The specific heat jump AC at the transition temperature is calculated from equation (p8|). It is shown in Figure 
||. Note that the jump is much larger than the rigid-lattice value of lA3jTp. The trend shown in Figure || can be 
explained by a rough argument correlating the sizes of AC/jTp and A(0)/&b?p. Simply put, if A(0)/kpTp is large 
then A(T) 2 will have a large slope at Tp. It has previously been noted experimentally^ that the order parameter has 
a BCS temperature dependence with A(0) and Tp treated as independent parameters. Some theoretical justification 
was recently provided for such a temperature dependence well away from TpEx Close to Tp the BCS form gives 

T \ 1/2 



A(T) ~ 1.74A(0) [ 1 - — ) (67) 



Within a BCS type of framework the specific heat discontinuity is given byfEl 



dA2 

AC ~~ P °lT 
Using A(0) = 1.76k B T RL and 7 = 2tt 2 Po /3 gives 



A(0) 2 , . 

= 3.03p o ^^. (68) 
1 P 1 p 



AC [T RL ^ 2 



1.43t2> V t p 



(69) 



This simple argument gives the correct trend that as the fluctuations increase the enhancement of the specific heat 
jump increases. 



D. Width of the one-dimensional critical region 

The width of the one-dimensional critical region is calculated from equation ( fiCf ) with the Ginzburg-Landau coef- 
ficients in the presence of the pseudogap. It is shown in Figure |j| normalized to the rigid-lattice value Ati£> = 0.8. 
The large reduction is very important because it means that even for weak interchain coupling, it may be possible 
for condition (Q) to be satisfied and for a mean-field treatment of a single chain functional, such as that used in this 
section, to be justified. 



VI. ESTIMATE OF THE PSEUDOGAP IN K0.3MOO3 



Optical conductivity, magnetic susceptibility and photoemission experiments all suggest that near Tp — 183 K 
there is a pseudogap in the density of states. 

ptical conductivity. Sadovskii has calculated the optical conductivity <j(u>) for the model introduced in ..Section 
For small frequencies <t(oj) is linear in to and has a peak at about uj ~ The data in ReferencesEjO then 
implies ip ~ 40 meV and ip/Tp ~ 2.5. On a less rigorous level ip can be estimated based on the analysis contained 
in the inset of Figure [|. If the single-chain mean-field transition temperature T < OATrl then ip ~ Tp^. Usuig the 
BCS relation (g) and the estimate A(0) ~ 80 meV for the zero-temperature gap from the optical conduct ivityc3 gives 
ip ~ 45 meV and ip/T P - 3. 
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Magnetic susceptibility. The data of ReferencesEn22 gives x(Tp)/x(^00K) ~ 0.5. Assuming that x(300K) ~ xo and 
using Figure |H gives tp/Tp ~ 2.4. 

Note that all of the above three estimates for tp/Tp are consistent with one another and are all in the regime where 
the fourth-order coefficent b is small (Figure ||) . n n 

Photoemission. Recent high resolution photoemssion measurementsc3L2l on K0.3M0O3 and (TaSe^I have several 
puzzling features: (1) There is a suppression of spectral weight over a large energy range (of the order of 200 meV for 
K0.3M0O3) near the Fermi energy. (2) The spectrum is very weakly temperature dependent. The suppression occurs 
even for T ~ 2Tp. (3) At Tp the spectrum does not just shift near Ep, due to the opening of the Peierls gap, but 
also at energies of order 0.5 eV from Ep. 

These features cannot be explained using the model presented in this paper. The photoemission data suggests 
that the pseudogap is about ip ~ 130 meV. Clearly this estimate is inconsistent with the estimates (tp ~ 40-50 
meV) given above from the optical conductivity and magnetic susceptibility. Furthermore, in the model presented 
here the pseudogap occurs only when £,\\(T) 3> vp/ ib, i.e.. fairly close to Tp. The temperature dependence of the 
Pauli spin susceptibility and the optical conductivityEaa suggest that the pseudogap disappears for T > 2Tp (in 
contrast to (2) above). Dardel et al.Ea speculate that the anomalous behaviour that they observe may arise because 
the photoemission intensity 1(E) might be related to the density of states p(E) by 1(E) = Zp(E) and thCj-auasi- 
particle weight Z vanishes due to Luttinger liquid effects. This suggestion has been examined critically by VoitEJ who 
concludes that the photoemission data is only quantitatively consistent with a Luttinger liquid picture if very strong 
long-range interactions are involved. Kopietz, Meden, and SchonhammerE3 have recently considered such models. 

VII. CONCLUSIONS 

In this paper a simple model has been used to illustrate some of the difficulties involved in constructing from 
microscopic theory a Ginzburg-Landau theory of the CDW transition. The main results are: (1) The large thermal 
lattice motion near the transition temperature produces a pseudogap in the density of .slates. (2) Perturbation theory 
diverges and gives unreliable results. This is illustrated by showing that a predictecO scaling relation between the 
specific heat and the temperature derivative of the susceptibility does not hold. (3) The pseudogap significantly alters 
the coefficients in the Ginzburg-Landau free energy. The result is that one-dimensional order parameter fluctuations 
are less important, making a mean-field treatment of the single-chain Ginzburg-Landau functional more reasonable. 

This work raises a number of questions and opportunities for future work, (a) The most important problem is that 
there is still no microscopic theory that can make reliable quantitative predictions about how dimensionless ratios 
such as A(0)/kpTp, AC/jTp, and £,q z Tp/vf depend on parameters such as Vp, the electron phonon coupling A, 
Tp and the interchain coupling, (b) Is the change of the sign of the fourth-order coefficient b of the single-chain 
Ginzburg-Landau functional for tp > l.lTp an important physical effect or merely a result of the simplicity of the 
model considered here? (c) Calculation of the contribution of ihp-sliding CDW to the optical conductivity in the 
presence of the short-range order associated with the pseudogapEJEj. 
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FIG. 3. Modification of the electronic specific heat C e (T) 
and the Pauli spin susceptibility x(T) by the pseudogap. Both 
are normalized to their values in the absence of the pseudogap. 
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FIG. 5. Dependence on the pseudogap of physical quan- 
tities associated with mean-field theory of a single chain. 
The plot shows the coherence length £o, the width of the 
one- dimensional critical region Atm, and the inverse of the 
specific heat jump AC. All quantities are normalized to their 
rigid-lattice values. For tp > 2.7T the coefficient b become; 
negative ijnd the transition will be first order (Section fv\ ) 
The large reduction of^57lly--below the rigid lattice value ol 
3ans tl at a mean-field treamtsnt of the single chain 
inzburg-La: idau functional may be justrfiejl 
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. 6lGhr The pseudogap duevtoYtherma^ lattice motion has^ 
(~\ a ^teni fcan t effect on the coefficients in the\Ginzburg-Landau 
free ensrgy (Q) for a single chain. \Tlre ratio of^the single-chain 
mean- ield transition temperature T\ and thq coefficients 
and c 1o their rigid lattice values (gWen\by (^|- O)) are shown 
as a function of the ratio of the pseudojmp tp to th^ tempera- 
ture. For tp > 2 .IT the coefficient b bkcon\es negativ^ and the 
transition will be first order (Section [eV|) . i\iset: Relationship 



(TtTq/Trl and iP/Trl- 



TABLE I. Comparison of experimental values for 
K0.3M0O3 of various dimensionless ratios with the predictions 
(rfitwo'rsimpls microscopic models. The three-dimensional 
trafMtkrnRLi iperature is Tp = 183 K. The zero-tempemture 
energjf gap A(0) is estimated from optical conductivity 
dataE3. A_Farni velocity of vf = 2 x 1ft 5 cm/sec was es- 
timated froiriband structure calculations!!!! AC is th« spe 
cific heat ju 1 up at the transition! 5 ] and jTp is the normal 
state" ejectroiic specific heat that has been calculated from 
the density^c fp*jtates estimated from magnetic susceptibility 
measurement 3Ej well above the transition temperature. I The 
longi^udiQd juheience lenglk £02 has been estimated n 
x-ray se^tte^ig experirnentsp. 0i /^pth models the dijien- 
Sionless rattag are independent pi any parameters, except for 
A(ffy/*kE.Tp in stehulz's mod^w^clyifwlescribed in Appendix 
[a|. The -rigid lattlc>vl^eor}My^r/olV(^. mean-field treat- 
ment of the 'single-chainSinzburg Landau functional (0) with 
the coefficients (P.-Jlll) 
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TABLE II. Parameters for several quasi-one dimensional 
materials. The observed transition temperature Tp is always 
much smaller than the rigid-lattice transition temperature 
Tb.l- The phonons near 2kp, which soften at the transition, 
can be treated classically since they have frequencies of the 
order of Q(0) (estimated from Raman and neutron scatter- 
ing) which is much smaller than Tp. The zero-temperature 
gap A(0), estimated from the peak in the optical absorption 
was used to calculate Trl {Trl = A(0)/1.76fcs). 





T P (K) 


A(0) (meV) 


Tp/T rl 


Q(0) (K) 


K0.3MOO3 


183 


90 a 


0.31 


80 b 


(TaSe 4 ) 2 I 


263 


200 c 


0.20 


130 d 


K 2 Pt(CN) 4 Br . 3 


120 e 


100 f 


0.18 


58 g 


TSeF-TCNQ 


29 


10 h 


0.42 





a Refill 

b J. P. Pouget, B. Hennion, C. Escribe-Filippini, and M. Sato, 
Phys.rAev. B 43, 8421 (1991) 



d S. Sugai, M.Sato, and S. Kurihara, Phys. Rev. B 32, 6809 



RefB 
S. Su 

(1985). 

c Complete ordering does not occurl 
f Ref|p 

SR □ 
h From activation energy of dc conductivity, Ref.tJ 
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APPENDIX A: SCHULZ'S MODEL 



For completeness an alternative microscopic model is discussed. SchulzE] considered only thermal fluctuations in 
the phase of the order parameter. He assumed that the temperature was sufficiently low that fluctuations in the 
amplitudG-Qf the order parameter were not significant. (However, in real materials the amplitude fluctuations are 
importantEil) . Fluctuations along the chain were treated exactly and the interchain interactions were treated in the 
mean- field approximation. He derived a free energy functional of the form (Jfy. The coefficients in a tetragonal crystal 
(a x = a y ) are 

2 2v F 

a = i-*p < A1 > 
B = ib (A2) 



C x = tj (A3) 



C = |fl ( A4 ) 



The coefficient A(T) vanishes at the three-dimensional transition temperature 

T 3D =(^) V2 . (A5) 
Thus the ratio A(0) /ksT^o is not a universal quantity. The specific heat jump at the transition is 



1 (A') 2 16ttT 3D 24 



where 7T = 2ttT/3vf<1x is the normal state electronic specific heat. The coherence length parallel to the chains is 



=X (A7) 



Equations (|A6|) and (A7) then give the ratios given in Table 



APPENDIX B: EVALUATION OF THE GINZBURG-LANDAU COEFFICIENTS 



The Ginzburg-Landau free energy functional ([!]) is related to the partition function Z by the functional integral 

Z 



Zn 



[#(*)] exp(-/3F[$) 



where Zq is the partition function in the absence of interactions. 

The phonon field in the Hamiltonian (To) is treated classically. For the Hamiltonian 



H = J dzy5>\z) 
the partition function is given by 



d 1 

Wircr 3 — + - (A(z)a+ + A(z)V_) 
oz 2 



*(*) 



A(z) 2 

XltVp 



Z = f dA(z)exp(-P f dz^^-) x 
J J Xttvf 



Texp 



dr / dz^(z, t) 



d 1 



(Bl) 



(B2) 



(B3) 
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The goal is to get this expression into a form-.comparable to (Bl) so the coefficients a, b, and c can be extracted. 

The linked cluster or cumulant expansionEZI can be used to rewrite the time-ordered product in flB3| ). In general if 
the Hamiltonian is separated according to 



H — Hq + Hi 



and 



(S) 



Texp ^- J* drffi(r)^ 



Z 
2o 



where 



)o denotes a thermal average with respect to Hq then the linked cluster theorem states that 

(S) = exp ((5) ,conn) - 1 



(B4) 



(B5) 



(B6) 



where (S)o,conn denotes the set of connected terms in the diag rammatic expansion of S. 
The important question is how to make the separation (B4)? 

A rigid-lattice treatment of the phonons neglects the effect of the thermal lattice motion on the electronic states. 
The Hamiltonian is separated according to 



Hi 



0,RL 



dz^(z) 



d 
oz 



(B7) 



Himl = \J dz^(z) (A(z)a+ + A(z)V_) *(*) 
The resulting free-energy functional is 

F RL [4>]= J dz {j^} + ^ U Te *P(-l j dZ ^ Z ) dT&(z,T)<T + V(z,T)-h.C.)J 



- 1 



(B8) 



(B9) 



,RL,conn 



where A(z) has been equated with 4>(zX-. in the functional integral. Expanding to fourth order in 4>{z) gives the 
rigid-lattice coefficients given by - |ll[)Ej. 

To improve on this rigid-lattice treatment we want to expand relat ive to a Hamiltonian which includes at least 
some of the effects of lattice fluctuations. The rigid-lattice expression (B9) is modified in the following way: 



F[4>] = J dz {j^) + ^ W Tex P(-^ / dzt{z)J* dT&(z,T)*+* 



(z,t) - h.c.) 



- 1 



where the expectation value of an operator A is defined by 

dvdv 



((A)) 



cxp(-vv*)Tv[Aexp(~(3H [v})] 



with the Hamiltonian Hq 



H [v] = J dzl&(z) 



iv F a 3 -^ + ^ip(va + + v*<j-) 



*(*) 



It is now possible to evaluate analytically the free energy functional (BIO). The Fourier transform 



is performed and the free energy (B1C) expanded to fourth order in (f>. The result is 
F[^]=^2,a{q)<t> q 4>* q + j b(q 1 ,q 2 ,q 3 )<f)* 1 <f> q2 <t> q3 (p qi + q2 - q3 



(BIO) 



(Bll) 



(B12) 



(B13) 



(B14) 
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where the coefficients a(q) and 6(91,92,93) are given by 

1 



a(q) 





f dk 


f dvdv* 









cr+G(k + 9, e„, v)a-G(k, e n ,v) 



(B15) 



6(91,92,93) = nT^2 



dvdv* 



dk - 

—Tr [a-G(k + qi,e n ,v)a-G(k + qi +q 2 ,e n ,v) 
a + G(k + 91 + 92 - 93, e n , v)a + G(k, e n ,v)} 



and the electronic Green's function G(k,e n ,v) is defined in equation ([mJ). 

Expanding a(q) in powers of 9 2 gives the Ginzburg-Landau coefficients a(T) and c(T): 



(B16) 



a(g) = a (T)+ c (T)9 2 + ... (B17) 
The fourth-order coefficient is b = 6(0, 0, 0). The coefficients a, 6, and c in are given in terms of ip in equations (|50|-^2|). 
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